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Abstract 

We present a detailed study of the dynamics of correlations in non-Mar kovian environments, ap- 
plying the hierarchy equations approach. This theoretical treatment is able to take the system-bath 
interaction into consideration carefully. It is shown that crosses and sudden changes of classical 
and quantum correlations can happen if we gradually reduce the strength of the interactions be- 
tween qubits. For some special initial states, sudden transitions between classical and quantum 
correlations even occur. 
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I. INTRODUCTION 



In the field of quantum information and computation, a great amount of meaningful in- 
vestigations have been done to explore the dynamics of two-level systems in environments 
with different properties jl-4|. Information extracted from a bipartite system comes to 
the center of our interests because in this simplest scenario we can not only conveniently 
introduce our theoretical calculations but also obtain the same results even with different 
definitions. In a two-qubit system, quantum mutual information has been defined as the 
largest amount of information which we can extract between the qubits j^]. This mutual 
information consists of two different parts, the classical correlation and the quantum correla- 
tion. Quantum correlation, which is equivalent to quantum discord in some special cases {(J, 
can also be divided as quantum entanglement and non-entanglement quantum correlation 
Q|. Many works have been focused on the different dynamics between quantum correlation 
and entanglement 

The most popular definitions of quantum information are discussed by Vedral, in which 
the total correlations (quantum mutual information) comes to the form as 8] 

Z(Pab) = S(p A ) + S(ps) - S(pab), (1) 

where pa{b) is the reduced density matrix of pab and S(p) represents the von Neumann 
entropy which can be calculated as S(p) = — SjAilog 2 Aj, in which A, is the eigenvalue of 
p. Classical correlation can be obtained from the measurements of subsystem B under the 
complete set of orthonormal projectors in the form 

n„ = sin 2 6>|1>(1| + cos 2 0|O)(O| + sin0cos0(e^|l)(O| + e~^|0> <1|) (2) 

and 

n ± = cos 2 6>|1)(1| + sin 2 6>|0>(0| - sin0cos0(e^|l)(O| + e _i *|0)(l|), (3) 
where |1) and |0) are the excited and ground states of the qubit. Then classical correlation 



can be expressed as 



C(pab) = max[S(p A ) - T lj q j S{p> A )] (4) 

0,<p 

with pP A = Tr B [(I(g)Il J )p ylB (I(8)n j )]/gj, qj = Tr^p^^I®!^)], and j =||, JL. Since the total 
correlations can also be expressed as 

Apab) = C{p AB ) + Q{pab) (5) 



then the quantum correlation is 



Q(p AB ) =I(p AB ) -C(pab). (6) 

With above basic equations, the left work is to give the values of the system matrix 
Pab- I n most scenarios we may find that the classical and quantum correlation would have 
mediocre behaviors, such as the result that classical correlation is always lager than quantum 
correlation. However, Many meaningful theoretical and experimental works have been done 
recently to demonstrate that with some special channels and a class of initial states, the 
classical correlation and quantum correlation will both remain constant in different regimes, 



10|4l2| 



and there is a sudden transition from the classical to the quantum decoherence 

In this work we demonstrate that in some common cases we can also obtain crosses 
of classical and quantum correlations and correlation sudden changes. By applying the 



hierarchy equations approach [J, [13] we set up the model of a two interacting qubits system 
in independent non-Markovian environments, and we also carefully take into account the 
presence of system-bath coherence, which has been proved to dramatically influence the 
dynamics of the quantum system. As the strength of the interaction between the two qubits 
decreases, we can find the crosses and sudden changes gradually appear. When there are 
not interactions at all, the changes will be the sharpest. 

II. THE MODEL AND THEORY 

We present a newly developed theoretical method called hierarchy equations approach 
to handle the system dynamics 0, [k^ . This method can conveniently treat the scenario in 
which there is a non-Markovian environment. We know that both the quantum system and 
the heat bath have characteristic frequencies. Most of theoretical treatments would suppose 
a fast bath to make the Born and ultrafast bath approximations valid {4]. However, in this 
hierarchy equations approach, a slow noise bath can be applied. This kind of environment 
will record part of the information from the system to the bath and then return some to the 
system. 

We set h = 1 throughout this report. The environment is divided into two independent 
baths interacting with two interacting qubits. Then we write the standard system Hamilton 
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as 



H s = e(a\ai + a\a 2 ) + ((a\ + ai)(a\ + a 2 ), (7) 

where e is the energy gap of the two-level system. The subscripts 1 and 2 represent the 
two qubits, a) and a are creation and annihilation operators. ( stands for the strength of 
the interactions between two qubits, so we can conveniently reduce its value until zero to 
observe the changes of the correlations. The heat bath consists of a set of oscillators, with 
the form 3, ^, 14] 



p 2 . 1 

Hb+sb = J^C^T + ~ S C ™3 ( a ™ + a m) x j, ( 8 ) 

j 3 m,j 

and Uj are coordinate, momentum, mass and frequency of the jth harmonic oscil- 
lator, respectively. Then Hb presents the Hamilton of the bath and Hsb is the interaction 
between system and environment. C m j stands for the coupling strength of the jth oscillator 
to the mth qubit. All of the system-bath coupling information can be found though the 
spectral density function J{oS) and in this work we model it as the Lorentzian cutoff 

7] represents the system-bath coupling strength, and 7 is the characteristic frequency of the 
heat bath. When the fast bath condition is fulfilled as ftry ^> e, we can choose the initial state 
in which the system and bath are independent. If fry ~ e it represents a non-Markovian 
environment and that is just our scenario. In the present work we choose e = 1.55 and 
7 = 45, and here 5 stands for an energy unit. With a Fourier transform of the spectral 
density function we have 

00 

F(t) = J^pfce-^W, (10) 

fc=0 

where jk is a Bosonic Matsubara frequency, with 70 = 7 and 7^ = 2nk/(3 when k > 1. The 
coefficients in front of the exponential function in Eq.[10]have values as Cq = ^(— i+cot ^7/2) 
and Cfc = 2?77o7fc//3 (7^ — 7^) when k > 1. The equation about the dynamics of the system 
using hierarchy equations approach, witch has the following form 4] 
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dp n {t) 



^ + ^^n mfc 7 fc )p n (t) - ^((o— -*^) 7 7-5Z— )[K», [Kn,pn( 
m=l k=0 m=l " k=0 ^ k 

2 K 

-i J2^2 Tlmk ^ CkVmPn ^k-^ t ) - C kPn mk -l(t)V m ) 
m=l fc=0 
2 K 

m=l fc=0 



where £ is the Liouvillian operator of the two- level system with Cp = [Hs,p], and V m = 
a^ m + a m with m — 1, 2. 

In Eq.[LTJ p n are called auxiliary density operators (ADOs) [3]. The subscript n is a 
multi- index which can be extended as n m ^, where m represents the number of qubits in 
the system and the k represents that of the Bosonic Matsubara frequencies. The notation 
n m fc±l refers to an increase and decrease of this index. Only when all the values of the multi- 
index are zero it stands for the physical system density operator. And if one of the notation 



n m k ± 1 has a value under zero or the sum of all the notations exceed a limit |13j] , the density 
operator is set to be zero. The values in the ADOs contain a deal of important information 
about the system-bath interaction, which has been proved to be able to dramatically change 
the dynamics of the system. With a carefully consideration of the system-bath coupling 
we find that it is a important factor which can greatly influence our results of the system 
evolution. Another advantage of this theoretical treatment is the convenience of applying 
the method into computer programs. It is quite easy to numerically calculate the dynamics 
of the system density matrix though a program. This hierarchy equations approach will be a 
useful treatment in scenarios when the interaction between the system and the environment 
is very strong. 

III. THE RESULT AND DISCUSSION 

In this work the environment is non-Markovian and much more common than the ones in 
Ref. . We set 7] = 0.35, 7 = 45 and j3 = 2.5/5. The initial state of the two-qubit system 
is a Bell state with \ip) = l/-\/2(|10) — |01)). As the system Hamilton only has diagonal and 
anti-diagonal elements, the system density matrix will evolve into X states [15], which can 
represent a great number of quantum systems. We use numerical method to calculate the 
mutual information, classical correlation and quantum correlation with enough precision. 
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FIG. 1: (color online) The green dashed, red dotted, blue dash-dotted and black solid lines represent 
the quantum correlations as functions of time with £ = 6, 0.75, 0.35 and zero, respectively. 
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FIG. 2: (color online) The black dotted, red solid and blue dash-dotted lines represent the mutual 
information, classical correlation and quantum correlation as functions of time, respectively. 

In order to show the correlation sudden changes, we set ( = 5, 0.75, 0.35 and respectively 
in Figured! Obviously, with the decrease of (, the quantum correlation becomes sharper and 
sharper. When ( = 0, a correlation sudden change occurs at certain time. For a complete 
view of this interesting change, we show the whole evolution of mutual information, classical 
correlation and quantum correlation with no interactions between the qubits in Fig.EJ The 
mutual information goes down monotonically as the system evolves. However, it is shown 
clearly that there are several crosses and sudden changes of the classical and quantum 
correlation. At first the quantum correlation shows a high speed of decay as a result of the 
interaction with the bath. Then the decay stops and the quantum correlation even shows a 
little increase, which clearly demonstrate the non-Markovian memory effects. However, the 
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classical decoherence speeds up and at around t = 1.7/5 these two correlations reach the 
same value. A cross immediately after the time shows the lager discord than the classical 
correlation during a short interval. Another cross emerges at t ~ 3.6/5. After this we 
may not find any crosses and the quantum correlation is then always larger than classical 
correlation. With careful observation we know that at both crosses and some other time 
points (such as t ~ 6.2/5 and t ~ 8.4/5) there are sudden changes of the correlations. 
These sudden changes catch interests because they may be results of some quantum phase 
transitions or the transition of classical and quantum decoherence 12J . We also find that the 
length of the intervals between these changes are all nearly 2.2/5. This phenomenon may 
indicate some mechanisms of the system evolution in a common non-Mar kovian environment. 

1.5, , 



• Eigenvalues 




o.o | —-^7 

1 2 3 4 5 

time(1/8) 

FIG. 3: (color online) The eigenvalues of the density operator pL after the maximization procedure 
as a function of time. The dotted line in the inserted figure at the top right corner is the derivative 
of the lower eigenvalue. 

In order to investigate the sudden changes, we calculate the eigenvalues of p 3 A after the 
maximization procedure in Eq.Hl Fig. [3] shows the first two sudden changes of the eigen- 
values. These two curves have exactly the opposite trends and both have two inflectional 
points, at which the derivative can become quite sharp, as we show in the inserted figure. 
Obviously the two points correspond to the sudden changes in the classical and quantum 
correlations. 

Next we concentrate on the influence of coupling strength. We choose ( = 0.55 and 
C = 0.15 for comparison which are shown in Fig.Hl The left figure shows a quite common 
dynamics of the correlations, with classical correlation always larger than quantum counter- 
part. When the strength of the interaction between the two qubits decreases to 0.15, crosses 
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FIG. 4: (color online) Mutual information, classical correlation and quantum correlation for dif- 
ferent coupling strength between the qubits. £ = 0.5<5 in the left figure and ( = 0.15 in the right 



one. 



begin to emerge as the quantum correlation rises and classical correlation goes down. How- 
ever, with careful observation we find that the changes are not so sharp as the ones in Fig.[T] 
where ( = 0. So we can infer that at the point of ( = some interesting phenomenon may 
happen such as the correlation sudden changes. 

Apart from the effects of the interactions between the two qubits, the initial states of 
the system also play an important role on the evolution of the correlations, as shown in 
Fig.0 We change the initial state from odd (|^) = l/v^dlO) — |01))) to even parity 
= 1/V2(|H) - |00))). The parameters are chosen to be the same as those in Fig. [3 
Sudden transitions between quantum and classical correlations can be found which has been 
demonstrated in other systems 111, [l2| . This kind of transition can occur several times in 



our system during the evolution such as at the time t = 1/5, t = 2/5 and so on. 
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FIG. 5: (color online) The black dotted, red solid and blue dash-dotted lines represent the mutual 
information, classical correlation and quantum correlation as functions of time, respectively, with 
an initial state of = 1/V2(|11) - |00». 

IV. CONCLUSION 

We study the evolution of the classical and quantum correlations of a two-qubit system 
in independent non-Markovian environments, applying the newly developed hierarchy equa- 
tions approach. The influence of initial state and coupling strength between the qubits are 
investigated. Crosses and sudden changes of classical and quantum correlations are illus- 
trated. We further show sudden transition between classical and quantum correlations for 
certain initial state. 
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